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Abstract 

By means of a recently introduced bihamiltonian structure for the homo- 
geneous Gaudin models, we find a new set of Separation Coordinates for 
the sl{r) case. 



1 Introduction 

In this paper we will discuss the Gaudin system with s/(r)-valued spins defined 
by the Hamiltonian 



Hg = J2Tt{A-A,), A^eslir) (1.1) 



i<j 

on the manifold M = sl{r)"' equipped with the standard product Lie-Poisson 
structure. We will refer to it, with a slight abuse of notation, as the homogeneous 
sl{r) XXX Gaudin system (as in P). 

The "conventional" approach to the integrability of this quite well studied 
problem is based on the Lax representation and the r-matrix theory (see, e.g., 
[IHIEH])- Fixing n distinct parameters a^, i = 1, . . . ,n one introduces the matrix 

1=1 

to be considered as an element of the Loop algebra s/(r)((A)). Along the Hamil- 
tonian flow deflned by (ll.lj) . the Lax matrix C{X) evolves according to a Lax 
equation 
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Thanks to the existence of an r-matrix for the Lax matrix ()1.2|) the spectral 
invariants 

jf") = Res| Tr(/:(A)"), i = l,...,n, a = 2,...,r, (1.3) 

are in involution. These integrals, together with the integrals of the motion 
associated with the invariance of the system under the global SL{r) action 
given by 

A, ^ GAG-\ (1.4) 

to be referred to as global gauge invariance, provide a complete set of constants 
of the motion for He- 

The separability of the Hamilton-Jacobi equations associated with the Gaudin 
Hamiltonian ()1.1|) was first studied [2H1 1211, 

kind of 

byproduct of the solution of the Bethe Ansatz equations associated with the 
quantum Gaudin system. Separability was then proved for the general case in 
|14| |2Zj . and (implicitly) framed within the theory of Algebraically Complete 
Integrable Systems in PlEj. In this scheme, it turns out that one can find a set 
of "algebro- geometrical" Darboux coordinates {d, Aj) as coordinates of a set of 
dn = r{r — l)(n — l)/2 distinguished points on spectral curve 

r(C,A)=Det(C-/:(A)), (1.5) 

(r — 1) 

whose genus is g = — - — ((n — 2)r + (r — 2)) . 

In the recent paper we have reconsidered the (homogeneous XXX) Gaudin 
model, and generalized to the case of an arbitrary Lie algebra g an alternative 
set of integrals of the motion for He (see, e.g. p]), introduced in the Hopf- 
algebraic approach to the integrability of the system. The distinguished feature 
of such integrals, which in the case of g = sl{2) are given by the very simple 
expressions 

ir, = Tr(^A,)', (1.6) 

1=1 

is that they are independent of the (fake) parameters entering the definition of 
the Lax matrix p.2p . 

These integrals were also introduced, in a different context JH], as Hamilto- 
nians of a notable class of Hamiltonian flows on the moduli space of + 3-sided 
polygons in M^, and later generalized in ^2] to the M.^ case^. This moduli space 
turns out to be a suitable Marsden-Weinstein Hamiltonian quotient of the Pois- 
son manifold sm(2)"+^ associated with the corresponding Gaudin model. The 
Hamiltonian flows associated with ()1.6|) were termed "bending flows" due to the 
following fact: if one draws, from a chosen vertex, the n possible diagonals of an 
n + 3-sided polygon, the flow associated with the Hamiltonian geometrically 
represents the bending of one side of the polygon along the k-th diagonal (the 
other side being kept fixed). 

-'^We thank J. Harnad for drawing our attention to these references. 
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The key point for the analysis performed in was the introduction, along 
with the standard Lie-Poisson structure P, of a particular second Poisson struc- 
ture, hereinafter called R. In the n = 3 case, this structure is defined by its 
Hamiltonian vector fields as follows: 





[^1 


< A2 = 




^3 = 


[^1 



dF dF 



dA2 


^ dA, 


dF 


dF 


dAi 


dA2 


dF 


dF 


dAi 


dA, 



9F dF ^ 



dA2 dAs ' dAs'' 
dF 

where are elements of sl{r) to be properly defined in Section |21 
dAi __ 
The Poisson pencil R — XP and the integrals p.fjj) fulfill standard Lenard- 

Magri relations, namely one can check that 

PdTr(A^) = 0, i = 1, . . . , n, RdTviAl) = 0, 

(1.8) 



RdTi 4^) = PdKa, a = 2, 

i=l 



For the general sl{r) case one can show that it is possible to find a sufficient 
number of polynomial functions in involution that provide a set of integrals of 
the motion alternative to the set defined by the Lax matrix ()1.2|) . They 
share with the integrals ()l.(j|l the property of being defined independently of the 
parameters entering the Lax matrix ()1.2|1 . 

In the last Section of [H] we addressed the problem of separability of such 
flows, in the framework of the so-called bihamiltonian approach to the SoV 
problem (see, e.g., [211 El I2DI El , and the references quoted therein). In partic- 
ular we solved it for the sl{2) case by means of explicit computations, showing 
that the separation coordinates associated with the pencil R — XP are ratio- 
nal functions of the natural coordinates {hi, Cj, fi) in s/(2)", and the separation 
relations are quadratic equations in these coordinates. 

In this paper we will solve the corresponding problem for the sl{r) case, 
with r arbitrary. This task will be accomplished by means of a careful mixing 
of techniques of the theory of Lax equations with r-matrix structure, and the 
theory of bihamiltonian systems such as those exposed in a series of papers by 
Gel'fand and Zakharevich |[15j and Magri and collaborators [SI 1201 • particular, 
we will make extensive use of (reflnements of) the results presented in |10|, ITT] 
concerning the Separation of Variables of systems with an arbitrary number of 
Lenard-Magri chains. The key points for the analysis we are going to develop 
in this paper are: 

1. It is possible to deform the Poisson tensor R into R in such a way that 

a) R is still compatible with P and restricts to the (generic) symplectic 
leaf S of P. 

b) The integrals K^"'^ deflned by the pencil R — XP are in involution also 
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w.r.t. R, although the recurrence relation they satisfy in relation with the 
new pencil R — XP are more complicated than the usual Lenard-Magri 
relations. 

2. It is possible to define n — 1 Lax matrices La linearly depending on a 
"spectral parameter" A such that any formal vector field X{X) which is 
"Hamiltonian w.r.t. the pencil R — XP induces a Lax equation on each of 
the matrices La- 

Thanks to the first property it will be possible to endow the generic symplectic 
leaf 5 of P with a special geometric structure, that is, a (1, 1) tensor with 
vanishing Nijenhuis torsion, whose "spectral" data will provide us with a set of 
separation coordinates for the H-J equations associated with Hq. 

Thanks to the second property, as well as other specific features of the de- 
formation R of R, (to be fully discussed in the core of the paper) we will be 
able to show that the separation relations are provided by the spectral curves 
of the matrices La- The distinguished feature of such a SoV scheme is that 
the separation coordinates are defined iteratively in subsets of dr = r(r — 1) 
coordinates, which are coordinates of a set of r(r — l)/2 points on a genus 

9 = {r-l){r-2)/2 

curve, irrespectively of the number n of "sites" of the Gaudin model. 

A word of warning: the set of coordinates defined in this way on S must 
be completed by a set of r(r — 1) coordinates associated with the global SL{r) 
invariance of the model, just like the set of integrals coming either from the Lax 
matrix ()1.2|) or from the construction discussed in [9^ must be supplemented 
by the set of integrals associated with the global gauge invariance of the model. 
However, since these integrals are associated with a sort of "cyclic" coordinates, 
they will trivially enter the H-J equations and the problem of separability. So 
in the core of the paper, we will often "forget" about them. 

The scheme of the paper is as follows: In Section |21 we will fix some conven- 
tions and notations to be used throughout the paper, review the results of 
to be used in the sequel, and introduce the Lax matrices La. In Section El we 
will briefiy recall the main points of the bihamiltonian scheme for SoV, and, 
in Section FH. II we will discuss how to apply such a picture to the sl{r) Gaudin 
models. Finally, in Section |3] we will give examples our constructions in the 
sl{2) and s/(3) case. In the last section we briefiy summarize the content of 
the paper and add a few comments. In order to simplify the presentation, we 
collected the proofs of some important but somewhat technical points in three 
Appendices. 
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2 The bihamiltonian structures and the Lax 
matrices 



Let Q be the Lie algebra sl{r). It is known that it (as well as any simple Lie 
algebra) can be identified with its dual, e.g., via the dual pairing given by the 
trace in the fundamental representation. In this paper we will constantly use 
such an identification. The Lie Poisson structure on M = is the one defined, 
in the natural coordinates {Ai, . . . , An} by its Hamiltonian vector fields: 

A=[A.,|£], (2.1) 

where, if X = (Xj,...,X„) represents a tangent vector to M, the elements 
dF 

-^-^ G s/(r), i = 1, . . . , n are those matrices defined by means of the expression 
of the Lie derivative of F w.r.t. X as 

" dF 
Liex(F) = 5^Tr(X,.-^). (2.2) 

1=1 * 

We will hereinafter denote the Poisson tensor associated with the Lie-Poisson 
natural bracket by P. From, e.g., [^12^]) we know that we can endow M with 
a multi-parameter family of Poisson structures which are compatible with the 
natural one (j2.ip . In [9j a further linear Poisson structure, to be denoted by R, 
has been introduced. It can be described as follows. 
We notice that relation ()2.1|) can be written as: 

Ai= 22 Pijk[Ak, -Q^], with Pijk = SikSij. (2.3) 
j,k=i ^ 

The new Poisson tensor R is analogously defined by the expression: 

dF 
dA. 

j,k=i ^ 
with "structure constants" given by 

rijk = {k — l)6ij5jk — 0{i-k)Sij + d{j-i)Sik + (^{i~j)^jk (2.5) 

where 6 is the usual Kronecker symbol and 6'(j) equals 1 if z > 0, and vanishes for 
i < 0. Explicitly, the Hamiltonian vector field associated by R with a function 
F is given by: 

4 r4 /• .N^^ \^ dF ^ dF dF^ 

= - 1)^- + 5: ^] + ^ - (2-6) 

* k=i+l " k=l * 



dF 



The following facts can be proven [5]: 
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Proposition 2.1 1. The pencil of bivectors R—XP is a bihamiltonian struc- 
ture on M , that is, R is a Poisson structure compatible with P. 

2. the functions = Tr{A'^~^^), a = l,,...,r — 1 are common Casimirs 
for R and P. The Lenard-Magri chains starting at 

H;^} = ^Tr{A:+'), a = 2,...,n 
a + 1 

provide us with further d = {n — l)r(r — l)/2 functionally independent 
integrals 

Hi^l, a = 2, . . . , n, p = 2, . . . , a + 1, a = 1, . . . , r - 1. (2.7) 

3. Taking into account the integrals associated with the global SL{r) in- 
variance of the model, that is the ring of functions generated by = 

■ Sr=i^«)' those integrals insure complete Liouville integrability of 
the model. 

□ 

Remarks. 

1. The Gaudin Hamiltonian is expressed in terms of the integrals ()2.7p 
as 

n n a—1 

a=2 a=2 6=1 

2. A convenient choice of the integrals associated with the global SL{r) in- 
variance can be done as follows. We pick the r — 1 independent elements 
Fhj^, . . . , Fh^_-^ associated with, say, the standard Cartan subalgebra of 
sl{r), and the Gel'fand-Cetlyn invariants, that is, the Casimirs of the 
nested subalgebras 

s/(2) C s/(3) C ■•• C s/(r), (2.8) 

under the map s/(r)" — > sl{r) sending the n-tuple {Ai, . . . , An} into the 
total sum, Atot = A^. 

For the sequel of the paper the following construction is crucial. Let us introduce 
n — 1 Lax matrices: 

a-1 

La = {\-{a-2))Aa + Y,^k a = 2,...,n, (2.9) 

k=l 

It holds: 
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Proposition 2.2 Let F be a smooth function on M and let us consider the 
pencil of vector fields 

= PxdF ■={R- \P)dF 

(we say that Xp is Hamiltonian w.r.t. the pencil P\). Then, along Xp, the 
matrices Li of eq. \2. .9)) evolve according to a Lax equation, 



with 



LiexdLa) = [La{\),Ma{\)] 



MJX) = (a-l-X)—+^ — 



(2.10) 



b=a+l 



Proof: Let us denote = i = 1, . . . ,n. The vector field Xp is explicitly 
given by: 



Lie^A(A,) = {P^dF)i = J^i^'ijk - Apyfc) [A 



k, 0(k\ 



^ ((/c - A - l)5ij5jk - 9{i - k)5ij + 0{i - i)Sik + 9{i - i)5jk) [Ak, a/c] 

i-l 

[Ak,ak - ttj] + 



k=\ 



Ai, {i-\ - l)ai + Qfc 

k=i+l 



Substituting in La{\) we get: 



'a-l 



Lie^f A [La] = {X - a + 2) i^[Ak, ak - aa] + 



,k=l 



a—1 I a—1 



j=\ \k=\ 

a-l 





N 


- A)aa + 


E 




k=a+l 


N 






)- 


k=j+l . 




N 




+ E 


k=a+l 





□ 

We can interpret this result by saying that we can associate with the homoge- 
neous n-particle Gaudin system a set of ra — 1 matrices depending on a parameter 
A, satisfying a Lax equation along the "formal" (i.e., depending on the param- 
eter A) flows of vector flelds that are Hamiltonian with respect to the pencil 
Px. 

Proposition 2.3 The coefficients K^\X) of the expansion in powers of n of 
the characteristic polynomial 

r-l 

deti^i - L,(A)) = fi' + 5^i^(:^)(A)/i^-"-i 



of every Lax matrix La{X) are polynomial Casimirs of the pencil P\ = R — 
XP. Moreover, along any vector field X associated with any of the non-trivial 
coefficient of such polynomial Casimir, all matrices La{X) evolve according to 
Lax equations 

Liex{La{X)) = [La{X),Ma{X)] 
for suitable matrices Ma{X). 

Proof. These assertion follow from the general theory of bihamiltonian pencils 
on loop algebras (see, e.g., [26^ and t25j). We sketch the proof for completeness, 
considering the equivalent set of spectral invariants Hm^ = l/(a + l)tr(Lm)°"'"^. 
To prove the first statement, we must show that, for any one-form f , that we 
can assume to be exact, v = dF we have 

{v,PxdH^J;^) = Q. 

Now, switching the action of the Poisson pencil on f = dF the LHS of this 
equation reads 

L^. = L^.{l/{a + l)Tr(L,(A))"+^) 

= ^ E Tr {l^{XY ■ L^. (L,(A)) ■ L„(A)"-^) (2.11) 

p=0 

= Tr {l: . L^.{lS) = tr (L^(A) • [L,(A), M^{X)]) = 

This proves the first assertion of the proposition, and, in particular, shows that 
all the vector fields xj^j, associated (say, via P) with the coefficients of the 
expansion 

i^^^ = E^S^' (2.12) 

are indeed bihamiltonian vector fields. 

To prove the second statement we notice, using a very simple trick well 
known to experts in the bihamiltonian theory of integrable system, that Xa], = 
PK^p_-^^ can be written as a Hamiltonian vector field w.r.t. the pencil, con- 
sidering the "truncated" polynomial (^X~PKa\X)^ , where (■)+ denotes the 

nonnegative part of the expansion in A. So we see that the bihamiltonian vector 
fields of the hierarchy are as well Hamiltonian vector fields w.r.t. the bihamil- 
tonian pencil Pa- The assertion then follows from Proposition 12.21 

□ 

Let us now focus our attention on a single Lax matrix, say L^; calling, for 
simplicity, Bg, = Ylb=i ^b, we have that the matrix 

= (A - a + 2)A-a + B-a 
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is a Lax matrix with spectral parameter that evolve according to Lax equations 
along the vector fields of the hierarchy. Clearly, the Poisson brackets induced on 
sl{r) by the map M M(2) defined by {Aa, Ba} are nothing but 
the Lie Poisson brackets on M^'^\ So, applying the formalism of IT] IHj.i.e.. 
according to the Sklyanin "magic recipe" jHOI, we can get, for every fixed a a 
set of canonical coordinates {^^, AJ^}. Actually, we shall do this in Section ITTl 

The point is that, to get a set of canonical coordinates for the whole systems, 
we have to compare the different sets of coordinates coming from the different 
Lax matrices La, a = 2, . . . ,n (and those coming from the global gauge invari- 
ance of the model). 

To solve this problem, we shall make use of the bihamiltonian structure of the 
problem, and, namely, frame the Gaudin systems within the so-called bihamil- 
tonian scheme for SoV. For the case of sl{2), we were able to solve the problem 
by means of straightforward computations. For the general case, we have to 
use some slightly more sophisticated ideas and techniques of the bihamiltonian 
theory, to be discussed in the next Section. 

3 Bihamiltonian geometry and Separation of 
Variables 

As we already remarked in the Introduction, a theory of Separation of variables 
based on the notions of bihamiltonian geometry has been quite recently intro- 
duced in the literature. The basic property of such a theoretical scheme which 
will enable us to solve the SoV problem of this paper can very simply stated as 
follows: 

Proposition 3.1 Let {M,Pi — XPq) be a bihamiltonian manifold and suppose 
that there exist functions f,g,Xf,\g, with A/ 7^ \g, (eventually, A/ and/or \g 
might be constant) satisfying 

Pidf = XfPodf, P,dg = \gPodg. (3.1) 

Then {f,g}, = {f,g}, = 0. 

Proof. The assertion easily follows from the equations 

{f,9}i = {df,Pidg) = Xg{f,g}o 
{9J}i = {dg,Pidf) = -Xf{f,g}o, 

□ 

In words, calling (with a slight abuse of language) a function / satisfying ()3.ip 
an eigenvector^^ of the pair Pi,Po relative to the eigenvalue!^ A/ as in |17j . 
this proposition simply says that eigenvectors belonging to different eigenspaces 
mutually commute. 

If the Poisson tensors Pi and Pq do not share the same image and kernel, 
then a complete set of eigenvectors cannot be found. This is a typical instance 
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in the Gel'fand-Zakharevich theory of bihamiltonian integrable systems, (and 
happens for the Poisson tensors R and P that we have considered so far). 

In general, the bihamihonian theory of SoV suggests to consider a suitable 
deformation Pi of Pi, such that it restricts to the (generic) symplectic leaves 
of P and it is still compatible with P. Upon restriction, the generic symplectic 
leaf S of P will be endowed with a regular bihamiltonian structure (that is, a 
bihamiltonian structure in which one element of the pencil is invertible). So, 
in the terminology of JT] the generic symplectic leaves of P are ujN manifolds, 
that is are symplectic manifolds (with symplectic form naturally induced by P) , 
endowed with a compatible Nijenhuis (or hereditary) tensor A^. In terms of the 
Poisson structures, the Nijenhuis tensor on the symplectic leaf S is defined by 

N = P,\s-Po%. 

To concoct out of Pi the suitable deformation Pi, one can adopt the following 
strategy: 

First one fixes a complete set Ci, . . . ,Ck of Casimirs of Pq, considers the first 
vector fields of the Lenard chains associated with Ca, i.e., 

Xa = PidCai a = 1, . . . , fc, 

and a a distribution Z, transversal to the symplectic leaves of Pq. For any basis 
Wi, . . . , Wk in Z, the matrix 

[Go\a,b = UewSCa) (3.2) 
is nonsingular (say on an open set U C M). So, the tensor defined by 

Pi = Pi - 5^ X, A [G^XtW, (3.3) 

a,b 

is well defined and restricts to the generic symplectic leaf S of P, since, by 
construction, PidCj = 0, j = 1, . . . , k. Notice that, if we define a new basis in 
Zhy 

Za = Y,[Go']tWt, so that LiezACb) = Sa,b, (3.4) 

b 

the expression of the deformed tensor Pi simplifies to 

Pl = Pl-J2XaAZa. (3.5) 
a 

We will call a basis of Z satisfying ()3.4p a normalized basis for the transversal 
distribution. 

The proof of the following Proposition can be found in fOl H] 

Proposition 3.2 Let (M, Pi — APq) be a 2n + k dimensional bihamiltonian 
manifold with corank{PQ) = k, and suppose that there exists an integrable dis- 
tribution Z C TM of dimension k, s.t.: 
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1. Z intersect transver sally the symplectic foliation of Pq. 

2. the space of functions invariant under Z is a Poisson subalgebra for the 
whole pencil (M, Pi — XPq)- 

Then, if Pi is the deformation of Pi defined by 4,'/. Pi — APq is still a Poisson 
pencil on M, and its restriction endows the generic symplectic leaves of Pq 
with the structure of a uN manifold. Furthermore, if Zi, . . . , Zj^ are a set of 
generators of Z, normalized w.r.t. a given complete set Ci, . . . ,Ck of Casimir 
functions of Po,condition\^ above translates into the equations: 

k 

Licz.Po = 0, LiczaPi = ^[Za, Xb]AZb, where Xa = PidCa, a = 1, . . . , fc. 

6=1 

(3.6) 
□ 

Definition 3.3 We say that a bihamiltonian manifold (M, Pi — APq), endowed 
with a transversal distribution Z satisfying the assumptions of Proposition \S.2 
admits an affine structure if it is possible to choose a complete set of Casimir of 
Pq, and a corresponding basis o/normalized flat generators Zi,, b = 1, . . . , k in Z 
such that, for every Casimir of the Poisson pencil H°'{X) and every b,c = 1, . . . k 
one has, in addition to Equation \3. 6|) 

Liez,UezXH\\)) = (3.7) 

The notion of affine structure for a biliamiltonian manifold was studied in ^T] 
in connection with the problem of the Stackel separability of a bihamiltonian 
system. For the purposes of the present paper, we remark that an affine Poisson 
pencil satisfies special properties, to be illustrated in the following. 

Let (M, Pi — APq) be a corank k affine bihamiltonian manifold, and let 
Za, a = 1, . . . , be a set of normalized fiat generators for the transversal distri- 
bution Z. Let us consider the polynomial Casimirs 

i7W(A) = A"-//o" + --- + K., 

and their deformations along the fiat generators, that is, the fc^ polynomials 

Dl{\) = Uez,H^^\\) = A""5," - Dl,X-"-' (3.8) 

The polynomials -D^(A) are invariant along Z, so that they can be considered as 
functions on the generic symplectic leaves of Pq. They are the building blocks 
of the bihamiltonian set-up for SoV for GZ systems. Indeed it holds: 

Proposition 3.4 Let A(A) be the determinant of the matrix of i3. then 

1. The roots Aj o/ A(A) satisfy 

PidXi = Aj Pq dXi, 
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2. Let -D(A) denote the classical adjoint matrix of D{X), and let [D{X)]a^c be 
non identically vanishing. Then any ratio p(A) := [D{X)]a,b/[D{X)]a,c of 
elements belonging to the a-th row of D{\), evaluated at the roots A, of 
A (A) satisfy the equation 

Pidp{\,) = KPdp{\i). (3.9) 
The proof of this Proposition is contained in Appendix A. 

□ 

3.1 Separation of Variables for the sl(r) Gaudin Systems 

In this subsection we will specialize the results of Section IHl and show how the 
bihamiltonian structure P\ = R — XP associated in Section |21 with the param- 
eter independent integrals of the Gaudin model provides a set of separation 
coordinates and relations for the H-J equations associated with Hq- 

The first step is to show that Px induces a uN manifold structure on the 
generic symplectic leaf S of P, that is, that the tensor R can be suitably de- 
formed. We consider in M = s/(r)" the n{r — 1) vector fields 

^"•=7^rT^' ^ = h---,n, « = l,...,r-l, (3.10) 

that is, the vector fields defined by their action on the n-tuple of matrices 
{Ai,...,An) by 

Lieiy^(Ai,...,A„) = (0,0,--- , e^,---,0), (3.11) 

i-th place 

where ea,r is the elementary matrix {ea,r)ij = Si,aSj,r- 

Proposition 3.5 The distribution Z spanned by the vector fields satisfies 
the hypotheses of Proposition 15*. M 

The proof of this Proposition is contained in Appendix B. 

□ 

To construct a set of fiat generators for Z, we can argue as follows. In the 
case of a single copy of sl{r), we normalize the with respect to the coefficients 
Ci, . . . , Cr-i of the characteristic polynomial of A. The normalization for the 
n site case is done site by site. Since the determinant of a matrix is a linear 
function of each of its entry, it is not difficult to realize that such normalized 
generators provide the GZ manifold M, R — XP with the structure of an 
affine GZ manifold, according to Definition 13.31 
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Let us now consider the Lax matrices 

a-l 

Lq = (A — a + 2)Aa + Ab, a = 2, . . . ,n, as well as Li = Ai. 

6=1 

Define Mj(A,^) = — Li{X), and denote their classical adjoint with Mj(A,^). 
The determinants of the matrices Mj define, thanks to Proposition I2.3| polyno- 
mial Casimirs Ka\X) for R — XP, via: 

r-l 

Det(M,(A,0) = r + 5^ir«(A)r~"-'. (3.12) 

a=l 

In particular, ii'i^'' are the common Casimirs of P and R, while i^'^")(A) are, for 
a = 2, . . . , n, the non trivial polynomial Casimirs of the pencil R — XP. 
Let us consider the matrices Dij defined by 

(A,); = Liez;(<(A)). (3.13) 

The proof of the following Proposition, which is based on a few elementary 
properties of the matrices (-Dij) is contained in Appendix C. 

Proposition 3.6 The determinant T>a{X) of the matrices [DaajjO' = 2, . . . ,n, 
factors as 

VaiX) = iX-a + 2y"'AaiX), (3.14) 

where Aa(A) is a monic polynomial of degree r(r — 1)/2. Let A*, a = 2, . . . ,n, s = 
1, . . . , r(r — l)/2, be the roots of Aa{X), and let us consider a row (say, the first) 
SaiX) = [DaaiX)) ^ ^ of the adjoint matrix Daa{X) of Daa{X). Finally, let he 
the functions obtained by evaluating in X = X^ the ratios 5r-2(A)/5r-i(A). Then, 
these [n — l)r(r — l)/2 pairs of functions {^l, A^} satisfy 

1) the Jacobi separation relations 

DetiMa{C,K)) = 0. 

2) the differential relations 

RdX: = XlPdXl RdC = KPdC- 

In particular, their brackets, ( say, with respect to the Lie Poisson structure P ), 
are of the separate form: 

{A^,e*}p = 5^*w:(c,A:), (3.15) 

where 9?a(Ca) A^) are functions of the two variables (^f, A^). 

□ 
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The meaning of this Proposition can be rephrased as follows. For every integer 
a = 2, . . . ,n we can construct a Lax matrix, whose characteristic polynomial 
gives us a family of Casimirs of the bihamiltonian pencil R — XP defined on 
the manifold M = s/(r)". Separated coordinates are constructed, according 
to the bihamiltonian scheme, by deforming such Casimirs along normalized 
generators of a suitable distribution Z defined in M. In particular, for each 
a = 2, . . . ,n, we can algebraically construct a "cluster" of (r(r — 1)) variables 
{-^a' ds=i,...r(r-i)/2 that are, in the terminology of [T7j, algebro-geometrical Ni- 
jenhuis coordinates, that is, satisfy properties 1) and 2) of Proposition 13.61 
To finish our job we have to: 

i) Discuss about the coordinates associated with the global gauge invariance of 
the Gaudin Systems 

ii) Explicitly construct, out of the coordinates found so far, a set of canonical 
separated coordinates (that is, a set of Darboux-Nijenhuis coordinates). 
Point i) can be solved as follows. One notices that any function ip depending 
only on the "global" matrix variable = Y17=i^i^ which is invariant along 
the distribution Z satisfies the differential relation 

Rdif = {n- l)Pdip. (3.16) 

In particular, this family includes the mutually commuting Hamiltonians of 
Gel'fand-Cetlyn type discussed in the Remark after Proposition 12.11 The prop- 
erty ()3.16|1 follows from the fact that they trivially satisfy the relation Rd(p = 
(n — l)Pd(f w.r.t. the undeformed pencil, and from the property that (f com- 
mutes with all the Hamiltonians of the hierarchy. Inside this ring of functions 
one can find a set of r(r — l)/2 canonical coordinates that complement the 
Gel'fand-Cetlyn Hamiltonians. Thanks to ()3.16|) they will have vanishing Pois- 
son brackets with the Nijenhuis coordinates of Proposition 13.61 

The solution to point ii) above can be simply done by means of a direct 
computation of the Poisson brackets between and A*. In particular, this 
computation will implicitly prove that these quantities are functionally inde- 
pendent. 

Proposition 3.7 The Poisson brackets, w.r.t. the Lie Poisson pencil P of the 
coordinates A^,^f defined above are given by 

{K,C} = iK'a + 2){K-a+l) (3.17) 

Proof. The proof of this proposition follows verbatim that of Theorem 1.3 in 
[T], to which we refer for full details. Indeed, the coordinates A^,^^ can be seen 
as common zeroes of the first row of the matrix Ma{X,^). So we can apply all 
the considerations of jl], the only difference being that the Poisson brackets of 
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the entries of the matrix Ma(A,^) are given by: 
{M:^(A,0,Mf(a,r^)} = 



tr 



a-l 



(A - a + 2)((7 - a + 2)Aa + {ekj^ii - eu^kj) 



r=l 



X-a 



[(A - a + l)(a - a + 2)(Mf (A, - M'^{\ + 



+ (A - a + 2){a -a + l){M'^{a,r^)5,k - Mi\a,7^)5u)\ 



(3.18) 



The presence of the factors (A — a + 2)(cr — a + 1) and (A — a + 2) (a — a + 1) is 
responsible for the factor (A^ — a + 2)(A^ — a + 1) in Eq. (j3.17p . 



□ 



4 Examples 

In this Section we will specialize the constructions presented in the paper for 
the cases of s/(2) and s/(3). 



4.1 The slf2) case 



Here we briefly reframe the explicit computations of the last Section of [H] 
within the formalism exposed in this paper. We consider the manifold M = 
s/(2)", endowed with the Poisson pencil R — XP of Section |21 It is explicitly 
parametrized in terms of the n matrices 



hi Ci 
ft -hi 



(4.1) 



The generic symplectic leaf of P is a 2n dimensional symplectic manifold, 
defined by the equations 



1 



and can be (generically) endowed with the 2n coordinates {hi, fi), i = I, 
A set of normalized transverse vector fields are given in this case by 

fi dci ' 



n. 



(4.2) 



(A - a + 2)ha + YTbJi hh (A - a + 2)ea + YTb=l ^6 



The matrices La are explicitly given by 
/ 

La{X) = 

\{X-a + 2) fa + EtZl fb -((A - a + 2)ha + Et=l h) J 



(4.3) 
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As canonical coordinates associated with the global SL{2) invariance one can 
choose the two functions 



45 V'l 



1 = 1 



The non trivial separation coordinates are gotten simply considering the zeroes 
Za of the elements [i^a]2,i, and the values /ia on these zeroes of the elements 
[-^a]2,2, normalized as in the previous Section. One sees that 



X]a:=1 fk 
fa 



+ (a — 2), a = 2, . . . ,n 



Shifting these values by the unessential term a — 2, we find that the separation 
coordinates are given, for a = 2, . . . , n, by 



fa 



^aha + X]fc=l 
XaiXa - 1) 



(4.4) 



They fulfill the separation relations 



2{XaiXa-l)) 



^ a—1 a—1 

-—TT^, {C'X + Tt{A^{J2 ^b))^^ + Tr((E ^^)')) • (4.5) 

" b=l 6=1 



In other words, the separation coordinates are coordinates of suitable points on 
the rational curves ()4.5|) . The corresponding Hamilton- Jacobi equations can be 
explicitly solved by means of algebraic functions. 



4.2 The sU3) case 



We consider the Poisson manifold M = s/(3)", endowed with the Poisson pencil 
R — XP and parametrized by the n matrices 



A,, 



\ 



% = 1, 



hi^i Ci^i 63^, 

y /3,i f2,i —h2,i I 

On this manifold the Poisson tensor P has 2n Casimirs: 

Q^ = ^Tr((A^)), Cf = iTr((A,)3) . = 1, 



n. 



n. 



(4.6) 



(4.7) 



The characteristic polynomials of the Lax matrices La = {X — a + 2)Aa + Y^l^-^ Ai, 
are expressed as 



r^{f,,x) = f,'-f,Ht\x)-H;t'W- 
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The transversal distribution Z is generated by the 2n flat generators: 
1 



d 
1 
d 



d d 

ifsAKi - ^2,i) + /2,i/l,i)^ h {f2,ihl,i - /3,iei,j)^ 

oe2,i oez,i 



d 



2,2 



de. 



d 



3,i 



3,i 



de 



2,2 



— /2,j/3,j(2/il,j — /;.2,j) + /2,j/l,i — fl,i^l,i 

The symplectic leaves of P are generically parametrized by matrices Ai of the 
form: 

^ hi,i ei,i $3,i ^ 

fl,i ^2,i ~ '^'2,i 

y /3,i /2,i — /i2,i y 

where $2,i and $3,i are suitable functions of the coordinates /ii^j, /i2,j, /2,i, /3,i, ei,i, 
parametrically depending on the Casimirs ()4.7|1 . 

The coordinates {A^,^^} can quite explicitly be found by means of the fol- 
lowing steps: ^ 

We consider the matrix Ma{X,i) = ^ — -^^^(A) and its adjoint Ma(A,^). We 
have to look for the common zeroes of the elements Ma(A, ^)3^i and Ma(A, ^)3^2, 
that is, for the common zeroes of 



Det 



^ -La{X)2,l e - ^a(A)2,2 ^ 



V 



-^a(A)3 (A)3,2 



Det 



J 



\ 



-L^X) 



3,1 



-LaiX) 



3,2 



(4.9) 

Taking into account the form of the vector fields and of the characteristic 
polynomial ()4.8|1 . we can identify the system ()4.9j) with 







CUez^M""^ + UeziH^^ = 



(4.10) 



where 



Det(M,(e,A)) = e'-ifJ")e-if(") 



(4.11) 

As we have noticed in Section ITT] , we can factor out (A — a + 2) from each line 
of this system, and consider, in matrix form the equivalent system: 



'-^2,2 '-^2,3 



'-^3,2 "-^3,2 



withG;^^^ = Lie^s-^^^VlA -a + 2). 



(4.12) 
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We notice that GJ^^ are polynomials in A of degree a — 1, so that the three 
zeroes A^, A^, A^ of the equation 



A„(A) = Det 



y '-^3,2 '-^3,2 j 



(4.13) 



are the compatibility condition for the system ()4.1()j) : the corresponding coor- 
dinates are thus given by, e.g., 

^ = -^^,3/^^,31.-.. (4.14) 

We remark that our procedure for finding separation coordinates exactly matches 
the one introduced, in the framework of r-matrix theory, in 0. 

Finally, defining (^^ = A^ — a + 2, a = n, . . . , n and considering the pairs 

{C, Pi} = C/CiC - 1), 3 = 1, 2, 3, a = 2, . . . , 

we see that the solution W of the (stationary) Hamilton- Jacobi equations of the 
s/(3) Gaudin model can be expressed as: 

where denotes the point (^^ , Ca ) ^^e genus g = I algebraic curve defined 
by equation ()4.11|) . and and denote, respectively, a suitable complete fam- 
ily of Gel'fand-Cetlyn Hamiltonians associated with the global SL{3) invariance 
of the model, and their canonically conjugated variables. 

Standard arguments show that the linearization of the flows associated with 
the mutually commuting Hamiltonians we have considered in this paper, and 
hence also the flow associated with the "physical" Hamiltonian of the s/(3) 
Gaudin system can be achieved by means of the Abel maps associated with the 
differentials (the last four being of the third kind) 

rf (c - i)r| ' (c - i)r^ ' C(C - i)q ' C(C - i)q ' ^ ' ' 



5 Conclusion and discussion 

In this paper we have reconsidered the SoV problem for the (homogeneous 
XXX) Gaudin systems based on the Lie algebras sl{r), from a particular stand- 
point. Namely, we considered the n-site system as a Gel'fand-Zakharevich sys- 
tem defined on the manifold M = s/(r)", with respect to the Poisson pencil 
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Px = R — XP, where P is the usual Lie Poisson bracket on M, while R, given 
by equations ()2.4|) and ()2.5p is a further linear Poisson structure on M that was 
introduced in [H]. 

We showed that the system admits a set of n — 1 "Lax" matrices, linear 
in the spectral parameter A that evolve according to Lax equations along any 
vector field that is Hamiltonian with respect to the Poisson pencil Pa- 

Thanks to this property, by using the bihamiltonian set-up for SoV, we 
managed to define a set of separating coordinates, quite explicitly given in 
equations ()4.4j) for the s/(2) case and by equations ()4.13|1 and ()4.14|1 for the 
s/(3) case. 

We notice that this set of coordinates provide an alternative set of separa- 
tion coordinates for the Hamilton- Jacobi equations associated with the Gaudin 
Hamiltonian w.r.t. the coordinates that can be found by means of the conven- 
tional approach based on the rational Lax matrix (|1.2j) . i,e, 

i=l 

In this respect, a few remarks are in order: 

First we notice that the existence of different sets of separation coordinates is 
to be ascribed to the super-integrability of the system. This is particularly clear 
in the bihamiltonian setting, where the separation coordinates are obtained, by 
means of the procedure outlined in Section El from the polynomial Casimirs 
of a Gel'fand-Zakharevich Poisson pencil. The general problem of the connec- 
tions between super-integrability and "multi-separability" is, to the best of our 
knowledge, still an open problem (see, ^] and the references quoted therein). 
In particular, the classification problem for these systems have been solved only 
for systems with a small number of degrees of freedom. The homogeneous XXX 
Gaudin models are systems with an arbitrarily high number of degrees of free- 
dom where the connection between super integrability and multi-separability 
happens, and their study might shed light on the structural properties of this 
phenomena. 

The second remark is the following. The SoV scheme based on the Lax 
matrix 15.11 leads to the definition of a divisor of degree 

dn = r{r~l){n-l)/2 

on the spectral curve i?(A,/i) = Det{fi — C{\)) (see, e.g., P^). It is not difficult 
to ascertain that the genus of this spectral curve is, for g = sl{r) 

^£ = ^^^((n-2)r + (r-2)), 

that is, it grows linearly with the number n of sites of the model. 

As we have shown in Section ITT] and exemplified in Section |3] for r = 2,3, 
the SoV scheme herewith outlined leads to consider the set of n — 1 Lax matrices 
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La{X) given by 



a-l 



Lai\) = i\-a + 2)Aa + Ba, Ba = Y,A,. 



(5.2) 



The separation coordinates parametrize sets of degree r(r — 1)/2 divisors on the 
spectral curves r'*(A,/i) = Det(/i — La(A)). We see that the genus of such curves 
is 



that is depends only on the rank r of the algebra, and not on the number n of 
sites, showing that the equations of motions can be explicitly solved by means 
of 9 functions of genus gY<=^ for all n's. 

This also implies the existence of canonical transformations between (suit- 
able open sets of) the degree dc Jacobian of the spectral curve -R(/i, A) associated 
with C{X) and the (corresponding open subsets) of the Cartesian product of the 
degree r(r — l)/2 Jacobians associated with the curves r"(A,/i)^. This simply 
follows from the fact that both the algebro-geometrical coordinates found from 
£(A) and those we discovered in this paper are canonical coordinates for the 
standard Lie-Poisson bracket Pq on sl{r)^, and, in particular, (together with 
those coordinates associated with the global SL{r) invariance) are Darboux 
coordinates for the restriction of Pq to its generic symplectic leaves. 

It is outside the size of this paper to fully discuss this issue here. However 
we think it is appropriate to display this transformation in the simple case of 
s/(2). This goes as follows. 

We recall that the n matrices Ai can be explicitly parametrized by means 
of 3n coordinates hi, Cj, fi, i = 1, . . . ,n: 



Let us denote with Ai,/ij the set of separation variables associated with the 
Lax matrix C{X) (see eq. fl5.H) ) of the n-particles s/(2)-Gaudin model and with 
Ci, Pi those associated with the bihamiltonian picture discussed in this paper. 

We can assume that the coordinates associated with the global SL{2) in- 
variance of the problem are the same; in particular, on the symplectic leaves of 
the Lie Poisson tensor Pq we have to consider the pair 



^7r. = (r-l)(r-2)/2. 



A, 



fi -hi 



with Lie-Poisson brackets given by 



{hi, ej}p = SijCj, {hi, fj}p 



5ij fj , {ci, fj}p — 2Sij hj 



i=i 

^This observation is due to B. Dubrovin. 




n 



Pi 
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According to the "Sklyanin" recipe, the A^, a = 2, . . . ,n are the zeroes of the 
(1, 2) entry of the rational Lax matrix ()1.2|) . i.e. the roots of the polynomiah 

n— 1 n— 1 n 

fe=0 ^j=l ^« fe=o j=l 

(5.3) 

Here the polynomials Sn-k-iidi, • • • , Oj, • • • , ctn) are the elementary symmetric 
polynomials in n — 1 letters bi, . . . , bn-i, defined by 

n— 1 n—l 

lli^- bj) = X;(-l)"''"'^n-,-i(&i, . . . , &n-i)A^ where so = 1, (5.4) 

j=l j=0 

evaluated for 6i = ai, . . . , bj_i = aj_i, bj = Oj+i, . . . , 6„_i = a„. 

We can express the "physical" coordinates fi in terms of the Aj as follows: 

Since the coordinates Q are rational functions of the fi alone, and namely, 

yi~] fk 

Ja 

we can explicitly find the transformation yielding the Q in terms of the Aj as: 

Y^( ^iA<^^-ai) e;:o<(-i)"-^'"Wi(A2,---,aj \ 

^ l^n,^,(a, - ai) ■ a^,(-l)-:'-is„_,_i(A2, . . . , A„) J ' 

for b = 2, . . . ,n, with, obviously 

Ci = Ai. (5.7) 

The variables fia, a = 2, . . . ,n are the values of the (1, 1) entry of the rational 
Lax matrix p.2|) for A = Aq, while the pa are be given by the values of the 
(1, 1) entry of the Lax matrix in A = Ca? divided by the normalizing factor 
Ca(Ca-l). Explicitly: 

a = 2,...,n, /ii = 

. , Aa - aj Ai 

(aha + E^J hj EJ=1 hj 

Pa = -rj-: rr a = 2, . . . , n, pi = 

The transformation of coordinates connecting the pi and the pi can be easily 
found noticing that they are connected to the coordinates hi by a linear trans- 
formation with coefficients depending on Q and Aj respectively. Consequently, 
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the transformation between the coordinates pi and hnear transformation 
with Aj depending coefficients: 

n 

Pi = ^ Aij(Ai, . . . , A„)/ij (5.8) 

i=i 

and it follows that it must be the lifting of the transformation defined by ( ()5.(ij) 
and ()5.7j) ) among the Q and the A^; 

n 

P. = Y.^{JT%^^o (5-9) 

where J is the Jacobian of the transformation ()5.6|) . 

Finally, we just mention some other problems which remain open. The first 
one is to compare our results with the picture of the generalized bending fiows of 

. The setting presented in that paper was aimed at providing a generalization 
of the previous paper [THj, and, in our setting, should be obtained by reduction 
to the submanifold of matrices Ai having rank equal to one. 

The second one is the application of the scheme herewith presented to the 
quantum sl{r) case. Preliminary results for r = 3 indicate that this should be a 
viable procedure for giving explicit expressions to the quantum integrals whose 
existence has been proven in [Hj. Work in both these directions is in progress. 
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Appendix A : Proof of Proposition 13.4 



The proof of Proposition 13.41 is divided into a couple of steps. We recall that we 
are considering a corank k affine bihamiltonian manifold (M, Pi — XPq), endowed 
with a transversal distribution Z, satisfying the requirements of Proposition l3.2l 
Za, a = 1, . . . , is a set of normalized fiat generators for Z, and 

i7W(A) = X^'^H^ + .-. + H^^, 

are polynomial Casimirs of Pi — XPq. Finally, the fc^ polynomials 

D,"(A) = Uez,H^^\X) = y^-5t - D^.X^^-' (A.l) 

are their deformations along the fiat generators. Finally, we recall the definition 

k 



a=l 
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Lemma 1 The actions of Pi and Pq on the deformation of the Casimirs of the 
pencil are related by the following formula: 

k 

PidDtiX) = \P,dDl{\) + Y,Dl{X)P,dDli. (A.2) 

c=l 

Proof. We limit ourselves to sketch the proof of this Proposition, which is 
essentially contained in Section 7 of [TT], although in a disguised form. We 
consider the characteristic property of a Casimir of the Poisson pencil, 

PxdH^''\X) = 

and derive it w.r.t. Zi,. We get: 

Uez,{Px) dH^''\X) + PxdDl{\) = 0. (A.3) 
Since Liez,(PA) = Y^AZb^^l] A with Xl = P^dH^ = PodH^, we see that 

[Zt,Xl] = Uez,{Xl) = UezSP^dHl) = -PodDl,. 
Thus eq. (IA.3|) takes the form 

PxdD-,{X) -Y^D^WPodDl, - Y,{[Zb,Xl],dH(^\X)) ■Z, = 0. (A.4) 

c c 

Let us consider the coefficient {[Zh, X^], dH''"-\\)) in the last sum. This is, by 
definition, 

Ue[z,,xi]iH^^\X)) = Uez^UexiiH^'^HX)) - Uex^Uez,iH^^\X)). 

Since {H^, H^°-\X)}o = only the second term is non identically vanishing, and 
equals -Uexi{D^{X)). 

Furthermore, thanks to the affinity of the GZ manifold, we see that all terms 
of the form Lie^^(iI'^(A)) identically vanish. 

So, we see that I A. 41 can be written as 

PxdD^iX) -Y,D:{X)PodDl, - Y,{Xl A Z,) ■ {dDl{X)) = 0, (A.5) 

c c 

which, in view of (|3.5p . yields the statement. 

□ 

Proposition 5.1 Let D be a k x k polynomial matrix of the form 

D«(A) = A-5» - X-^-'Dl, Dl^^,a,b = l,...,k, (A.6) 

where the ^ are smooth independent functions on a bihamiltonian manifold 
M, satisfying equation HA.S^) . Then: 
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1. Its determinant A(A) has the form 

A{X) = - AiX"-' + ■ ■ ■ + A„ (A.7) 
with V = J2a and satisfies 

PidA{X) = XPodA{X) + A{X)PodAi. (A.8) 

2. The roots Xi of A{X) satisfy 

PidXi = XidXi. (A.9) 

3. Let -D(A) denote the classical adjoint matrix of D{X), and let [D{X)]a,c be 
non identically vanishing. Then any ratio p(A) := [D{X)]a^h/[D{X)]a^c of 
elements belonging to the a-th row of D{X), evaluated at the roots Xi of 
A (A) satisfy the equation 

Pidp{X,) = \Podp{Xi). (A.IO) 



Proof. The power expansion ()A.6|1 simply states that the (a, a) entry of 
is a monic degree Ua polynomial, while all other entries in the a-th row are of 
degree not exceeding Ua — 1- We preliminarily notice that 

k k 

A(A) =1[d:{X)+0{X'^-'), whence A, = Y.Dl,. (A.ll) 

a=l a=l 



We multiply the matrix equation ()A.2|) say, on the left, by the classical adjoint 
D, to get 

Y,Pi{D:d{X)Dl{X)) -Y,^PoD{X)ldDl{X) = 

' (A.12) 

j2m>^)m>^)PodDi,. 

c,d 

Recalling that 'Ec^cW^'dW = ^adA{X) and 

J^DcWdD^iX) = TiD{X)dD = dA{X), 

a,c 

taking the trace of the matrix equation lA. 121 and taking into account (jA.llj) we 
get the proof of the first item. 

To prove item 7^ 2, we first notice that, for any function / on M, the 
evaluation of a polynomial (or rational) function F{X) in the parameter A, whose 
coefficients are themselves functions on M gives rise to a new function F{f) on 
M. Its differential can be written as follows: 



, dFjX) 
\=f dX 



df 
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To clarify the notations, the first term in the RHS of the above equations means 
the differential of F{X), with A taken as a parameter, then evaluated for A = /, 
and the second the partial derivative of F{X) w.r.t. A, subsequently evaluated 
for A = /. Keeping this proviso in mind, we consider now F[X) = A(A), and 
f = Xi, i = 1, . . . ,n. We have: 

dA{X) 



d(A(A,)) = rfA(A) 



+ 



A=A,: 



dX 



dXi , 



(A. 13) 



where (iA(A) = — -^^ ^d/S.j. Taking into account the relation ()A.13|) . we get 



= (Pi - XiPo)dA{Xi) = A{Xi)PodAi 
which implies the assertion, since 

A(Ai) = 0, while 



c'A(A) 



dX 
dA{X) 



A=A, 



iP,-\Po)dK, (A.14) 



dX 



A=Ai 



7^0, 



thanks to the fact that being the coefficients Aa^j functionally independent, the 
roots are generically simple. 

The proof of the third assertion is basically contained in the proof of Propo- 
sition 8.4 of ^I]. We limit ourselves to sketch it. 

By using the relations ()A.2|) and ()A.8|) . together with the defining relation 
i5°(A)-D^(A) = SadA{X), one arrives at the matrix equation 



[PidD: - >^PodD^^)Dl = A{X)Po{dA,6^ - dDl,). 



(A.15) 



If cr denotes one row of the adjoint matrix -D(A), we can rewrite the above 
equation as ^ 

(Picr - XPo(t)D{X) = A(A) ■ X 

where X is the corresponding row of the RHS of ()A.15|) . If we consider the 
normalized row p = cr/aj, we see that , since 

pD(A) = A(A)o-/a„ 

it holds 

(Pip-APop)^ = A(A)-F, (A.16) 
for some suitable Y whose form is irrelevant here. Evaluating this equation for 



A = Aj, we see that 



(Pip-APop)D(A) 







A=A, 



(A.17) 



Taking into account that A has simple eigenvalues, we see that each row 
{Pidp — XPodp) , evaluated at A = Aj must be proportional to the corresponding 
row of -D(A), that is, there must exist vector fields X' such that 



{P,dp-XPodp)\ 



x=x. 



X'.p 



A=Ai- 



Since one element of p is normalized to 1, we thus see that X' must vanish, 
whence the thesis. 
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□ 



Appendix B : Proof of Proposition 13.5 



The key point is the following observation on the (ordinary) Lie-Poisson brackets 

on a single copy of M = sl{r). The Poisson bracket of two functions F,G on 

OF dG dF dG 

M, is given by {F, G} = Tr(— ■ [A, — ]) = -Tr(A ■ [— , — ]). Let denote 

the ij-th entry of A, and consider the family of r — 1 vector fields on M defined 
d 

by = — - — , a = l,...,r — 1, as well as the distribution Z C TM defined 
by the W. 

We notice that differentials of functions vanishing along Z admit a very 

simple matrix representation. Indeed W"' is represented via its action on the 

matrix A as the elementary matrix ea,r, having 1 in the a-th place of the last 

fdF\ ' dF 
column. So Liew^iF) = iff — — = 0, a = 1, . . . , r — 1, i.e., iff — — lies in 

the lower maximal parabolic subalgebra p_ of sl{r). 

Let now W denote any element in Z, and let F, G be functions such that 
LiezF = lAezG = 0, and let us compute Liez{{F,G}). Thanks to the Leibniz 
property of the Lie derivative and the fact that Z is a constant vector field we 
have that 

/ [f)F f)G' 

Uew{{F, G}) = -Tr Uew{A) 



which vanishes as well since p_ is indeed a Lie subalgebra of sl{r). 

In the case of the n-particle sl{r) Gaudin model, whose phase space is 
parametrized by n matrices Ai, we consider the family of n ■ (r — 1) vector 
fields defined by 

Uewt{Ai,...,An) = {0,0,-- - , e^,---,0). (B.2) 

i-th place 

The distribution Z generated by these vector fields is generically transversal to 
the symplectic leaves of the Lie-Poisson product structure on sl{r)^ . We now 
prove that the space of functions vanishing along Z is a Poisson subalgebra for 
any affine Poisson tensor Q. The brackets {F, G}q = {dF, QdG) are given by 
the multiple sum 

N / r,r^ / N 



{F.G} 



Q 




Ak, 



dG 
dA~ 



dG 
'dA, 



where (Jk denote constant matrices. Noticing that the differentials of functions F 

, dF dF dF 

vanishing along Z are represented by n-tuples of matrices dF 



'dAidA2 dAr, 
, A^, we see that the Lie derivatives Liew/ajF, G}q are 



. , dF 

with — — ep_, 1 = 1 

dAi 

given by multiple sums of terms like those of Eq. ()B.1|) , and so vanish whenever 
Lieiyc-F = lAew^G = 0. 
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□ 



Appendix C : Proof of Proposition 13.6 



The proof of Proposition IH.fjl follows from a few elementary but important facts 
following from the definitions of the Casimirs of the Poisson pencil K^'^\\) and 
of the normalized transversal vector fields Z", a = l,...,r — 1, i = 1, . . . ,n. 
Recall that we defined matrices Dij by 



where the polynomial Casimirs Ka"^ are defined by 

r-l 



(C.l^ 



Det(M,(A,0) = r + E^'-^(^)^''"""'- 



(C.2) 



a=l 



One has: 



1. The matrix (-Dii) is the identity. This trivially follows form our choice of 
Li = A^. 

2. The n{r — 1) x n{r — 1) matrix of the deformations of the Casimirs w.r.t. 
the transversal vector fields has the following block form: 



Dn ■ ■ • 

D21 D22 



k3 









(C.3) 



This follows form the fact that, for j > i and every a, Lie^QLj(A) = 
0. Thanks to point 1 above, we now consider the non trivial matrices 
0) a = 2, . . . , n, and the corresponding fields Z^. 

Taking into account that: 

i) Liez^.(Det(M,(A,0)) = -Tr(M„(A, OLiez,"(^a)); 

ii) Liez^iLa) = (A - a + 2)Liez«(v4a); 

iii) The determinants of the diagonal blocks DetDaa = 'D"'{X) are monic 
polynomials of degree r(r + l)/2 — 1 in A; 

We can factorize T''^(A) as 

V'iX) = {\-a + 2Y~'A''{\) (C.4) 

where A"(A) is a monic polynomial of degree r(r — l)/2. 



27 



3. Thanks to the lower diagonal block form of the matrix D of equation ()C.3|) . 
every diagonal block Daa satisfies Proposition (^. So its determinant 
satisfy, according to Proposition O (-R - XP)dV{X) = V{X)PdV'l, and 
thanks to the factorization property ()(y.4jl . we have 

{R - XP)dA%X) = A"(A)PrfA^. (C.5) 



We now recall that the Casimir functions of P are given by the highest 
order terms of the expansion of the Casimirs of the Poisson pencil Ka\X) 
in powers of A. If we call 

(C.6) 



the matrix of the deformations of the Casimirs of P with respect to the cor- 
responding vector fields introduced in eq. ()3.10p . noticing that Lie^y^i^'^'^Q 
vanishes for j ^ i, we see that the normalized generators and the "constant" 
ones are related by 



n. 



Thus, considering only the nontrivial indexes a = 2, 



LieH/,.Det(M„(e, A)) = G^'^Liez.Det(M,(e, A)), 



(C.7) 



we can argue as follows. Since G"'^ are independent of A, we see that the 
common solutions (^a, A^) of the two sets of r — 1-tuple of equations 



and 



LieH/2Det(M,(e,A)) 

LieH/„.Det(M,(e,A)) 
Liez2Det(M,(e,A)) 

^ LiezjDet(M„(e,A)) 

coincide for every (fixed) a = 2, . . . , n. 

Expanding the RHS of the equations of the system ()C.9|) as 



r-l 



r-1 



Liez„.Det(M,(e,A)) = r-^-^Liez^.if^) = ^^^-^(Z},)^, 

(3=1 (3=1 



(C.9) 



(C.IO) 



we see that: 
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a) The roots of X>"(A), introduced in eq. ()C.4|) are those values of A for which 

the r — 1 equations, defined for a = 2, . . . , n, 

j2c-'-"m'i = o (c.ii) 

admit solutions. In particular, the roots of A'^(A) define non trivial ele- 
ments A*, s = 1, . . . , r(r — l)/2. 

b) The values Q corresponding to the roots A^ of A" (A) are given by suitably 

normalized elements of the adjoint matrix Da, evaluated at A = A^. 

Thus, from the bihamiltonian theory, we can conclude that the only non van- 
ishing Poisson brackets between such functions admit the separate form: 

{xiC}p = 5^*w:(c, a:), {k,C}~r = 6''6a,Kv:{CK) (c.12) 

We now consider eq. ()(I8jl . taking into account the observation that the 
pairs (Ag,^*) are solutions of this system as well. We notice that from the 
definition of W^, Lievi/^Det(Ma(^, A)) is nothing but the determinant of the 
minor of Ma(^, A) relative to the a,r entry. Since the r(r — l)/2 pairs (^f. A"), 
for every fixed a, annihilate the r — 1 minors of the matrix Ma{^, A) relative to the 
entries (1, r), . . . (r — 1, r), they annihilate the minor relative to the (r, r) entry 
as well. Hence they annihilate the last row of the adjoint matrix of Ma{^,X), 
and so satisfy the characteristic equation 

Det(M,(e:,A:)) =0, a = 2,...,n, s = 1, . . . , r(r - l)/2. 
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